On the comparison of volumes of quantum states 
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Abstract 

This paper aims to study the a- volume of /C, an arbitrary subset of the set of 
N X N density matrices. The a-volume is a generahzation of the Hilbert-Schmidt 
volume and the volume induced by partial trace. We obtain two-side estimates for 
the a-volume of IC in terms of its Hilbert-Schmidt volume. The analogous estimates 
between the Bures volume and the a-volume are also established. We employ our 
results to obtain bounds for the a-volume of the sets of separable quantum states 
and of states with positive partial transpose (PPT). Hence, our asymptotic results 
provide answers for questions listed on page 9 in [1] for large N in the sense of 
a-volume. 

PACS numbers: 02.40.Ft, 03.65.Db, 03.65.Ud, 03.67.Mn 

1 Introduction 

Recent development of quantum information and quantum computation theory has 
attracted considerable attention on the geometry of quantum states, which studies the 
geometry of the set and/or subsets of quantum states on the (composite) Hilbert space T-t 
(of complex dimension N). Among those subsets, the set of all quantum states ("D), the 
set of separable quantum states (5), the set of entangled quantum states {S := V\ S), 
and the set of quantum states with positive partial transpose iVVT) are of particular 
interest. In applications (see [2, 3]), entangled quantum states play fundamental roles 
because of the Einstein- Podolsky- Rosen (EPR) correlations [4] (see also [5]). Hence it 
is important to study (1) the probability of finding separable quantum states within 
(2) the necessary and/or sufficient conditions of a quantum state being separable 
(or entangled). Unfortunately, it turns out that the second problem is hard [6]. One 
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important and powerful necessary condition for separability is the well-known Peres- 
Horodecki positive partial transpose (PPT) criterion [7]. However, it is known that the 
Peres- Horodecki PPT criterion is not sufficient in general (except in very special cases; 
for instance, [8, 9, 10]), and entangled quantum states with PPT have been constructed 
[11]. Thus, a natural question is (3) how precise is the Peres-Horodecki PPT criterion 
as tools to detect the separability? Another question regarding to the Peres-Horodecki 
PPT criterion is (4) is the Peres-Horodecki PPT criterion precise as a tool to detect 
entanglement? To answer questions (1), (3) and (4), one often needs to study the size of 
S, S, VV1~, and V (for various relevant measures of size). 

In literature, measures on V have one common feature: unitary invariance. These 
measures can be written as the product of measures on simplex (of eigenvalues) and 
the measure on the manifold (of eigenvectors). Important measures on T> include, for 
instance, the Hilbert-Schmidt measure (V^s), the Bures measure (Vb), and the measure 
induced by partial trace on composite systems. The Hilbert-Schmidt measure is induced 
by the Hilbert-Schmidt metric which induces the flat and Euclidean geometry into V. 
Hence, a lot of known techniques, such as, techniques from geometric functional analysis 
and convex geometry, can be used to estimate the Hilbert-Schmidt volume of (convex) 
subsets of V, e.g., S and WT in [12, 13]. These two papers provided answers to 
questions (1) and (3) for large in the sense of Hilbert-Schmidt volume. In [14], the 
author obtained similar results in the sense of Bures volume, by comparing the Bures 
volume with the Hilbert-Schmidt volume. 

The present paper strives to answer questions (1) and (3) for large in the sense of 
the a- volume (see section 2 for its definition). The a- volume is a natural generalization 
of the Hilbert-Schmidt volume and the volume induced by partial trace on composite 
systems. The latter one has appeared in many places and attracted a lot of attention 
[15, 16, 17, 18, 19, 20, 21]. It is inspired by the important procedure of purification 
and the (up to a multiplicative constant) unique measure on the space of pure states 
of the (complex) Hilbert space "H. Any quantum state p on "H may be obtained by 
partial tracing of some pure quantum state over the auxiliary subsystem "H' (with 

dimension K > N), where |0) is a unit (column) vector in the composite Hilbert space 
"H ® "H'. The space of pure states of "H is isomorphic with the complex projective space 
CP^~^, and hence the unique measure on it, Pn{-), is the one induced by the Haar 
measure on the unitary group U{N). Therefore, to obtain Vn,k, the volume induced by 
partial trace on the composite Hilbert space H ^ H', one chooses the natural measure 
Pnk{-) on the space of pure states of "H 0^', and considers its push-forward induced by 
the the operation of partial trace [21]. When K = N, V^^k coincides with the Hilbert- 
Schmidt measure. Thanks to the work of Sommers and Zyczkowski [20, 21], one knows 
the precise mathematical formula of Vn,k and some statistical properties of Vn,k- In 
particular, they calculated the exact value of VN,Ki1^) for all A^ and K. However, the 



2 



calculation of the exact values of Vn,k{<S) and Vn,k{VVT) seems to be difficult because 
the geometry of these sets is not very well understood. These quantities can be used to 
measure the probabilities of separability and of PPT within V. In the present paper, we 
provide estimates of Vn,k{,S) and Vn,k{'PVT), as well as the a- volumes of S and WT . 
Our results show that the probability of finding separable states in T> is extremely small 
and so the Peres-Horodecki PPT criterion is not precise for (even moderate) large N, in 
the sense of a-volume. 

The present paper is organized as follows. In section 2, we introduce some necessary 
mathematical background and notations; all the results cited there are known or 
straightforward. Our main results are presented in section 3, that is, we compare the 
a-volume with the Hilbert-Schmidt volume and the Bures volume. Section 4 contains 
estimates of the a-volume of S and WT, and these new estimates supply solutions for 
questions (1) and (3) for large N. Moreover, we provide some numerical examples to 
show the effectiveness of our new estimates. 

2 Mathematical Background and Notations 

We work on the (complex) Hilbert space "H = C"^^ (g) C^^ ■ ■ ■ ® C^" with n > 2 and 
Di > 2 for all z = 1, 2, ■ ■ ■ , n. The (complex) dimension of equals N = D1D2 ■ ■ ■ Dn- 
Any quantum state on % can be represented as a density matrix, i.e., a N x N positive 
(semi) definite matrix with trace 1. Here the trace of a matrix is the sum of its diagonal 
elements. We use V = ViTi) to denote the set of all quantum states on "H. The set of 
separable quantum states on Ti [22] is denoted by S, that is, 

S = S(n) := conv{pi ® ■ ■ ■ ® p„. Pi G V{C^')}. 

The complement of S within V is the set of entangled quantum states, i.e., S := V \ S . 
Note that the entangled quantum states play crucial roles in quantum information 
and quantum computations. Both V and S are convex subsets with (real) dimension 
d = N^-l. 

For any quantum state p G there are some unitary matrix U G U{N) and some 
diagonal matrix A = diag(Ai, ■ ■ ■ , Xn) with (Ai, ■ ■ ■ , \n) G A, such that, p = UAW . 
Hereafter, A'^ is the complex conjugate of the matrix A, U{N) denotes the A^-th unitary 
group, that is, the set of all A^ x A^ unitary matrices. A refers to the regular simplex 

in M^, i.e., A = |(Ai, ■ ■ ■ , A^v) G : A^ > 0, ^JI^ A^ = l| . Obviously, the eigenvalue 

decomposition p = UMJ'^ is not unique. One can change the order in which the 
eigenvalues of p occur by using unitary permutation matrices. Without loss of generality, 
one can choose (Ai, ■ ■ ■ ,\n) G Ai, the subset of A with the order Ai > ■ ■ ■ > A^v > 0. 
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This corresponds to divide A into A^! parts, and to pick just one of them. The reason for 
estabhshing an order in the eigenvalues of p is to avoid duphcation in calculating volume 
of subsets of V. In the case of non-degenerate spectrum, i.e., Ai > A2 > ■ ■ ■ > \n, it is 
easy to verify that 

p = f/Af/t = UBkB^U\ 

where A = diag{\i, ■ ■ ■ , Atv) and B = diag{zi, ■ ■ ■ , zn) with \zj\ = 1 for j = 1, ■ ■ ■ , TV. 
That is, the matrix U is determined up to the N arbitrary phases entering B, and the orbit 
will be the coset space J-"^ = W(A^)/[W(1)]^. Again to avoid duplication in calculating 
volume of subsets of V, we will use measure on J-"^ instead of measure on hl{N). Note 
that if degeneracies occur in the spectrum of p, the matrix B need not to be diagonal in 
order to commute with A, however, the degenerate spectrum can be ignored because its 
total measure is 0. 

Measures on "D, such as, the Hilbert-Schmidt measure, the Bures measure, and the 
measure induced by partial trace, are invariant under conjugation by a unitary matrix. 
Moreover, they all have the product form: dv x ^7, where v are some measures on the 
simplex A and 7 is the invariant measure on J-"^. The 7 measure may be written as 

^7 = W 2Re{U-^dU)ijIm{U-^dU)ij, 

l<i<j<N 

where dU is the variation of U such that U,U + dU G U{N). Note that the 7 measure 
is the unique measure (up to a multiplicative constant) induced by Haar measure on the 
unitary group IA{N). It is known that 7(J-'^), the total 7 measure of J-"^ (see [23]), is 
equal to 

ln^\N{N-l')n 

tN 



where E(Ar) = n;=ir(j). (1) 



E{N) 

Here T{x) = t^~^e~^ dt is the Gamma function. We point out that 

-^<r(x)<-, or ^<r(l + a;)<l, for all x G (0, 1), (2) 

VX X v 

where ^ » 1.12917 [24]. 

For all a > 0, we define the a- volume, Vq,, as 

TV 

dVe. = n n - ^i)' 

i=l l<i<j<TV 

Here, for simplicity, we let 

TV TV 

dA = 5o(5^A,-l)n^A„ 
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where 60 is the dirac measure at 0. The a- volume of V can be calculated as follows (see 
[21, 23, 25, 26]) 

The Hilbert-Schmidt distance between any two states p, a G is defined as 



Dnsip, = \\P -ctWhs = ct)2). 

This (natural) metric equips V with the flat, Euclidean geometry on V, but it is not 
monotone [27]. It induces the Hilbert-Schmidt measure Vhs^ which equals to \/NVi 
[21, 23], namely, 

dVHs = v^ n (A..-A,)'rfArf7- 

l<i<j<N 

Note that the Hilbert-Schmidt measure is same as the usual (translation invariant) 
Lebesgue measure on V. Moreover, the Hilbert-Schmidt measure does not have 
singularities. By formula (4), the precise value of the Hilbert-Schmidt volume of V 
[23] equals 

V„s(V) = (2n)^ Vn (5) 



In later sections, we are interested in the ratio ^^^^j for IC C V, which is equal to the 



ratio 1^ j . Hence, it is convenient to ignore the constant VN for our analysis, and we 



Vhs{IC) 

-.UlVJllO, VVC Cll C llllj^i ^3Dlj^U. Ill Ull^ 1 dljlVJ 

keep using Vhs instead of Vi to represent the 1-measure, i.e., 

dVHS= n i^^-^J?dAd^. (6) 

l<i<j<N 

Another special case of the a- volume is the volume induced by partial tracing of the 
composite quantum system l-L^l-C. Here H. and l-i' are and K dimensional Hilbert 
spaces respectively. Without loss of generality, we assume K > N . The partial tracing 
over l-i' gives a reduced density matrix of size N x N. Any state p on "H (g> "H' can be 
expressed uniquely as 

N K 

i,3 a, 13 

where {ci}^^ and {fa}a=i ^^e canonical bases of "H and Ti' respectively. Define the 
partial trace over "H' as 

K 

= Treip), where p^- = ^pi^j^ for i,j = !,■■■ ,N. 

(5=1 



The measure induced by partial trace is an alternative way to derive measures on 
D. In fact, the partial trace process allows us to view states on "H as a (pure) state on 
(much higher) dimensional space "H ® T-i' ■ Then the measures induced by partial trace 
may be considered as a projection of the [NK — 1) dimensional simplex of eigenvalues 
into simplex of [N — 1) dimension A [21]. It takes the form 

N 

dVN,K = n n - ^^•)' 

i=l l<i<j<N 

where (Ai, • • ■ , Xn) G Ai. This is the a-volume with a = K — N+1, a natural number 
bigger than or equal to 1, and thus, (integer) a-measure is an induced measure on V. In 
particular, Vn^n is just the Hilbert-Schmidt measure. In other words, the Hilbert-Schmidt 
measure can be viewed as an induced measure on V. 

The Bures measure Vb(-) [16, 25, 28] can be formulated as follows: 

2-iV-iV^ 

V Ai ■ ■ ■ Xn i^t<J<N + 
where (Ai, ■ ■ ■ , X^) G Ai. It is induced by the Bures metric -D_b(-, ■) [29, 30], which takes 



the following form: for any two states p,a ^ V, Db{p,o') = \jl — Itr ^f^Jpo'^. The 
Bures metric is proven to be equivalent to [31] 
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DB(p,a) = sup f v yt^^w^ - \f^^¥^f\ 

where the supremum runs over all the projection measurements {IPi}^i, i.e., all P, are 
rank 1 projections satisfying ^^2!i=\ — IdN- Employing the projection measurement 
{^i}iLi to a state p results in the z-th outcome with probability tr(Pjp), z = 1, ■ ■ ■ , A^. 

Note that ( Xlili [V^^i^iP) ^ ^) ^ the Hellinger distance, which measures the 

statistical distinguishability between the discrete probability distributions {tr{Fip)}^^ 
and {tr(PjCr)}^^. Hence, the Bures distance provides the measurement of statistical 
distinguishability between p and a. We point out that the Bures metric is Riemannian 
but not flat. Moreover, it is monotone, that is, it does not increase under the quantum 
channels (completely positive, trace preserving maps) [32]. The precise value of the Bures 
volume of V was obtained in [28], 



A. X ^-V\l---\N,^t}<N + r(iVV2) ' 
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which is the rf- dimensional volume of the (i-dimensional hemisphere with radius |. (There 
is no satisfactory explanation for this mysterious fact.) Noticed that the Bures volume is 
not translation invariant and has singularities on the boundary of V, where at least one 
of the Aj's are 0. ( If two or more of them are 0, then some denominators in (7) vanish.) 
As in the case of Hilbert-Schmidt volume, it is convenient to ignore the constant term 

2 2 in (IVb- Hereafter, we keep using (IVb to represent the measure 

VAi---Ajv A, + Aj 

The primary goal in [14] is to compare VRb(/C,P) with VRh5'(/C,P) for any subset 
K. in V. Here the Bures volume radii ratio, VRb(/C,£), for two subsets IC,C gV, was 
defined as [14] 



vR.(K,£) = i^; 

where d = N"^ — 1 and the Hilbert-Schmidt volume radii ratio was defined as 

These quantities aim to compare the volume ratio of /C to C Similarly, we define here 
the a-volume radii ratio, VRq,(/C,£), as 

This can be used as a measure of the relative size of /C to £ in the a-volume sense, 
and does have geometric meanings. This paper strives to estimate VRq,(/C, V) from both 
upper and below in terms of VR//5(/C, V). Our proofs rely on the Stirling approximation, 
which can be written as 

We refer the readers to the references [14, 16, 21, 23, 25, 28, 33] for more detailed 
background and for motivation. 

3 a-volume VS Hilbert-Schmidt and Bures volume 



This section aims to compare the a-volume of subsets of quantum states in terms of 
their Hilbert-Schmidt volume and Bures volume. Hereafter, we let JC be any (Borel) 
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measurable subset of ViTi), the set of quantum states on (complex) Hilbert space "H. 
The following lemma, which is independent of possible tensor product structures of H, 
is our main tool to estimate the a-volume in terms of the Hilbert-Schmidt volume. The 
same approach will give corresponding comparison results for real Hilbert space "H. 

Lemma 1 Let a > be a (fixed) constant independent of the dimension ofH. Let be 
any (complex) Hilbert space with (complex) dimension N and JC be any measurable subset 
mV{H). 

(i): For all p > a > 1 , 

VnsiJCr (v^^^iV)^ ' < K(/C) < iV(^-")^ VnsifC), (9) 
(a): For all < p < a < 1, 

N('-»)^ Vhs{IC) < Vo^ilC) < VHsilCy (V)^ (10) 

Remark. From formula (4), for > 0, one has 

1 — ■» / / rV — r)\ TV 7- TV 7- / TV 7- -|\ \ \ 



T{{^)N + NiN-l))fi V Vl-P 



The above formula, as a function of p, can be extended to < 0; however there are 
singularities whose exact locations depend on and a. 

Proof. We only prove the case of a > 1. Similar approach will lead to the argument for 
the case of a < 1. 

To that end, let /i : A ^ M be /i(Ai,--- , Xn) = n^Ii^""^- Note that (9A, the 
boundary of the simplex A, consists of sequences for which one or more of the Aj's equal 
to 0. Therefore, /i(Ai, ■ ■ ■ ,Xn) is always on dA if a > 1, and is strictly positive in 
the interior of the simplex A. On the other hand, /i(Ai, ■ ■ ■ ,Xn) has unique critical 
point (1/A^, ■ ■ ■ , 1/A^) in the interior of simplex A by the Lagrange multiplier method. 
The compactness of h{Xi,--- , Aat) implies that h{Xi,--- ,X]sf) must have a maximum 
inside the interior of the simplex A for a > 1, and hence the (unique) critical point 
(1/A^, ■ ■ ■ , 1/iV) must be the (only) maximizer of h{Xi, ■ ■ ■ , A^r) on A. Therefore, 

N 

hiX,, . . . , A^) = n < iV(i-")^, for a>l. (11) 

i=l 
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'1<*<J<A^ 

By substituting g into the second integral of (14), one has, for p > a > 1 

N X 1-p 



(13) 



By formula (3), the a- volume of /C can be calculated as follows: 

N 

Va{IC)= l[xr' n (A.-A,)2 ciAci7. (12) 

'^'^ i=l l<i<j<N 

Therefore, formulas (6) and (11) imply 

Va{)C) < iV(i-")^y^5(/C), for a>l. 

Now let us prove the lower bound for a > 1. By (12) 

VailC) = [l[xr' n i>^^-XJ?dAd^= f f g dAd^, 

■^^ i=l l<i<j<N 

where, to reduce the clutter, we denoted 

TV 

^?(Ai,--- ,A;v)=n^r' n i^^-^.f"'"' 

1=1 l<i<j<N 

/(Ai,--- ,atv)= n i^^-^jf'- 

l<i<j<Ar 

For any p > 1, we employ the Holder inequality (see [34]) to (13) and get 

V^{]C)>(^jj-pdAdj^ (^j^g^hdAdj^ . (14) 
By formula (6), and substituting / into the first integral of (14), one gets 



(15) 



/ n iA.-A,f n^^^'''^^^^ 

i<j<i<Af 1=1 ' 

>( [ n \X.-X,\'flxF^"^dAdS~' ={V^{V)y-', (16) 

^ l<i<j<N 1=1 ^ ^ 

where the inequality follows K d T) and p > a > 1. The lower bound is then an 
immediate consequence of inequalities (15) and (16). 

Remark. More generally, one can compare V"q(/C) in terms of Va(/C), for all > 0. 
The above approach gives the following corollary. 
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Corollary 1 Let a,P > be (fixed) constants independent of the dimension ofH. Let 
"H be any (complex) Hilbert space with (complex) dimension N and K, be any measurable 
subset in ViTi). 

(i): For all p > ^ > 1, 

(a): For all < p < ^ < 1, 

The following lemma is simple but important for the proof of our main results. 
Lemma 2 Let a > be a constant independent of N. Then 



Remark. A direct consequence of this lemma is: let a, /3 > be two constants 
independent of N, then 



Proof. Now by formula (4), one gets 

VhsCD) r{N^-N + aN) / r(j) 



For each j = 1, ■ ■ ■ , A^, and by a > 1, one has 

m _ ix2x...(j-i) ^ 1 



r(a + j - 1) r(a) X a X (a + 1) X ■ ■ ■ (a + j - 2) ~ r(a) ' 
Combining with Stirling approximation formula (8), one has, for a > 1 



1/{N^-1) 



lim ^ii^i-^ < 1. 



(^^) ) 
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Now we bound 



1/d 

from below. By induction, for all j = 1, ■ ■ ■ , A^, and by 




T{j + a-l) - r(A^ + a-l)' 



and hence 




(18) 



Together with Stirling approximation formula (8) and formula (17), one 



lim 

N^oo 



( 




>1, 



which completes the proof. 



Theorem 1 Let a > be a (fixed) constant independent of the dimension ofH. There 
are universal constants (independent of dimension N, hut depending on a) Ci,c'i > 0, 
such that for any Hilhert space H and any subset K, gV, 



Moreover, the bounds are optimal in general. 

Remark. In fact, as — >■ oo, Ci{N,a) — 1 (and c[{N,a) — )■ 1). For small dimension 
N, one can precisely calculate Ci{N,a) (and c'i{N,a)) by formula (19), if a is given. In 
fact, one can let 



Moreover, by Lemma 1, Ci{N,a) > 1 (and Ci{N,a) < 1) for all a > and all A^. 
Let us consider n = (C^)®" with = 2". Then, if a = 2, Ci(4,2) ^ 1.19861, 
Ci(8,2) ^ 1.10785, Ci(16,2) ^ 1.05706, Ci(32,2) ^ 1.02958 which seems to decrease 
to 1 as A^ increases; while if a = 1/2, c;(4, 1/2) ^ 0.873608, c;(8,l/2) 0.936271, 
4(16, 1/2) ^ 0.968159, 0^(32, 1/2) ^ 0.984148 which seems to increase to 1 if A^ increases. 
Similar phenomenon happens to other fixed a, and this strongly suggests that Ci(4, a) 
for a > 1 or c'i{A, a) for a < 1 should work for all 2" > 4 (and, in fact, all other integers 
A^ > 4). 

Proof. We only prove the case of a > 1. Now dividing Va(T>) from both sides of the 
upper bound of inequality (9), then 



(t): for a>l, VR„(/C,P) < Ci VRhs{}C,V), 
(ii): /or < a < 1, VR«(/C,I?) > c'l VRh5(/C,I?). 




VhsOQ 
Vhs{V) 



) 



( 



) 
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Taking d-th root from both sides, then 



VR„(/C, V) < VRnsiJC, V) V^^^^) ' ^^^^ 

By Lemma 2, for a > 1, 

i'-JL 1, j 

and hence, there is a universal constant Ci > such that VRa(/C,P) < CiVRh5(/C, ^^)- 
To see the optimahty of the upper bound, we let /Ci = tl^ + (1 — t)pmax for < t < 1, 

i.e., 

]Ct=\uXU^:X = dmg{^j^ + tAi, ■ ■ ■ , ^ + ,(Ai, ■ ■ ■ A^) G A and ?7 G W(iV)| . 
Let Ztv be as in (1). 

Clearly VRH5(/Ct,P) = t by homogeneity of the Hilbert-Schmidt volume. Now we 
estimate Va{K^t) from below. In fact 



\4(/Ci) = z^t^'-' / n + n - ^^)" n ^^^^ 



1 i=l ^ ^ l<i<j<N 3=1 

As ^ + tXj > tXj for all j = 1, ■ ■ ■ , A^, and a > 1, one has 

„ N N-1 

K(^*)>t"^'^^+"^^"^^)^A^/ n\"^' n (A.-A,)^n^A,=t^^-^+^("-i)K.(p). 

Equivalently, 

N(a-l) 

Hence, if t > exp{—c'{{N'^ — l)/N) for some constant c" > 0, then 

VRa{]Ct,V) > exp{-c[{a - l))t = exp{-c[{a - l))VRiy5(/Ct, P). 

Part (i) guarantees that YRaiJ^t^T^) < CiYRusiJ^t)'^)- So the upper bound can be 
obtained for a > 1, and it is optimal in general. 
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Theorem 2 Let a > be a (fixed) constant independent of the dimension ofH. There 
are universal constants (independent of N, but depending on a) Ci,C[ > 0, such that for 
any Hilbert space % and any subset K, (I'D (let Q = VRhs'(/C, "C) ), 



(i) : for a>l, VR^(/C,P) > Ci VR^5(/C,P) exp 

(ii) : /or < a < 1, YRa{^,V) < C[ VR^silC^V) exp 



(1-Q)lnln(e/C) 



(l-a)lnlii(e/C) 
Af2-1 



Proof. We only prove (i). By the lower bound of (9), for all p > a > 1, 
Dividing Va{T>) from both sides, one has for all p > a > 1, 



Vaiv) - \Vhs{v)) V v;(P) ; V v^iv) 

Equivalently, by taking the d-th root from both sides, 

VR.(/C,P) > (VR„.,(/C.D))'(^) [-^VJvr) ■ 

We first bound I TTT-tin — I from below. To that end, we let r = t{N, Q = 



= Nir!(^e/() ^ function depending on and ( = WRusiJ^iT^)- This implies that 

p{N, Q = y5^, and p{N, Q ^ a a.s N oo. Moreover, < tN = < 1, and hence 
Tj < 1 for aU J = 1, ■ ■ ■ , iV - 1. 

By formula (4), one has 



VriV) _ T{N^-N + aN) fT{j + T-r 



V^{V) T{N^ -N + tN) jJ^\T{j + a-r 
Formula (2) and induction show that for all j = 2, ■ ■ ■ , A^, 

r(j + r-i) ^ r(r + l) ^ 1 



r(j + a - 1) - r(a + 1) - r(a + 1) ' 
because < r < 1. Thus, again by (2) 



n Vr(j + a-l)J - \T{a + 1)J T{a) " r [1(1 + a)] 
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N ■ 



This further imphes, hj p = > 1 and r 



1 



7Vln(e/C) ' 

hm n f > hm — — 



■J 



^ V f a .^^ /(I -a)lnln(e/C) 

/(I -a)lnln(e/C)\ , , 

hmexp ^ 1 . M - (21) 



Stirhng approximation formula (8) imphes that 



hm — — = 1. 



N^oo 



rViN^ - N + aN)\ 
\T{m - N + tN) ) 

Together with inequahty (21), one has 



(Ym^> exp ((1^4^^ I . (22) 



Inequahty (9) imphes that 



hm ^^/^ > hm A^^^^ > hm exp ^ ' ^ = 1. (23) 

Since T\n( = jv(i+in^i/c)) — W ^"^^ < r < 1, one has 

hm VRhs(/C,P)p-" = hm exp( rlnC ) > hm exp( i^—^ ] = 1. (24) 



Combining (20), (22), (23), and (24), one gets VR„(/C, V)>ciVR%s{'C, V) exp 
for some universal constant Ci > 0. 

Remark. A slightly more precise calculation shows that 



(l-a)liiln(e/C) 



Moreover, from the above proof, one gets for any fixed a > 1, the constant ci(A^, a) — )■ 1 
as — oo. Numerical results strongly suggest that, for any fixed a > 1, ci{N,a) 
increases to 1 as A^ increases, and ci(4, a) should work for all A^ > 4. For instance, if 
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we consider H = (C^)®" with = 2", then our proof yields that ci(4,2) ^ 0.219056, 
ci(8,2) ^ 0.374655, Ci(16,2) ^ 0.539382, Ci(32,2) ^ 0.686509, Ci(64,2) ^ 0.800189, 
ci(128,2) ^ 0.878831, ci(256,2) ^ 0.929124, and ci(512,2) ^ 0.959595. Similar 
phenomenon happens to the case of < a < 1, namely C'i{N,a) decrease to 1 as 
increases. For example, if a = 0.5, then C;(4,0.5) ^ 1.43603, C((8,0.5) ^ 1.24359, 
C;(16,0.5) ^ 1.13648, C;(32,0.5) ^ 1.07664, (64, 0.5) ^ 1.04292, and hence, C[{4:,a) 
should work for all > 4. 

As a consequence of Corollary 1, one can prove the following result, whose proof is 
similar to Theorems 1 and 2. 

Corollary 2 Let a,f3>0 be (fixed) constants independent of the dimension of T-i. Let 
% he any (complex) Hilbert space with (complex) dimension N and K, he any measurable 
subset in V {T-i). Let C, =YB.i3{lC,'D), Omax = maxjl, anc? = inin{l, |}. There 
exist universal constants (independent of N, but depending on a and /3) 02,6*2 >0, such 
that, 

Remark. As in Theorem 1, ]Ct attains the upper bound ii a > (3, and attains the lower 
bound if a < p. Therefore, for a > f3 the upper bound is optimal in general, and for 
a < P the lower bound is optimal in general. Theorems 1 and 2 are special cases of this 
corollary with /3 = 1. 

Theorems 1 and 2 in [14] compare the Bures volume of /C in terms of the Hilbert- 
Schmidt volume. The following theorem compares the Bures volume with the a-volume. 

Theorem 3 Let a > he a (fixed) constant independent of the dimension of Ti. Let 
"H be any (complex) Hilbert space with (complex) dimension N and JC be any measurable 
subset in Vil-L). There exist universal constants (independent of N, but depending on a) 
03,6*3 > 0, such that 

where ^ = VR„(/C,I?). 



15 



Proof. Similar to the proof of Lemma 1, rii<i<i<Af x +x ^'ttains the minimum at 
Ai = 1/iV for alH = 1, ■ ■ ■ , A^. Hence, 



> - 



l<i<j<Af 

Vi(/C). 



Recall d = N"^ — 1. Dividing Vb{T>) from both sides and taking the d-ih root, one gets 
that for some universal constant c" > 0, 

VRb(K.D) > l-j VR.(K.D)(^^ 

> c"VRi(/C,I?). (25) 
where the second inequality follows Lemma 2 and the limit [14] 

The lower bound is an immediate consequence of Corollary 2 and inequality (25). The 
upper bound is an immediate consequence of Theorems 1, 2 and Theorem 2 in [14]. 

Remark. If a > i, the lower bound is optimal in general, and if a < 1, the upper bound 
is optimal in general. If a G [1/2, 1], then there exist universal constants C3, C3 > 0, such 
that 

c'3VR«(/C,P) < VRb(/C,P) < [VR«(/C,I))]i/2exp i^^^^^j^ 



and the bounds are optimal in general. Moreover, there exist universal constants 

^3' "-^3 



Co, Co > 0, such that 



4'VRi(/C,P) < VR5(/C,P) < C'i [VRH5(/C,I?)]^/^exp i^^^^ 

with C = VRi^s(/C,P) and the bounds are optimal in general. This inequality improves 
results in [14]. 
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In fact, it was proved that ct < VRsilCt, V) <Ctiit<'^ for some universal constants 
c,C > [14]. Recall Kt = tV + {1 — t)pmax- Therefore, Kt attains the lower bound for 
a > 1/2 by a similar calculation in the proof of Theorem 1. 

To see the optimality of the upper bound for a < 1, we let < t < 1, and 

IC' = {UXU^ : X =diag(l -t + tAi,tA2, - ■■ ,tA^), {Xu---Xn) G Ai and f/ G W(iV)}. 

Recall Ai is the chamber of A with order Ai > ■■■ > Ajy. For e'^^"^"^^ < t < |, 



one has VRB(/C*,r') > VRi^5(/C*, P) where c, c > are two constants [14]. To 
prove the optimality of the upper bound for a < 1, it is enough to prove that 
VRhs{}C\V) > cVRa{JC\V) for some constant c > 0. 

To that end, we first bound Vq,(/C*) from above. By formula (3), one has 

N 



■^^ij=2 2<i<j<N 
N N~l 

x(l + tAi - t)""^n(^ -t + tXi- tX,f W dX 



r 

i=2 j=i 



As 1 - t < 1 + tAi - t < 1, < 1 - t + tAi - tA,- < 1, and a < 1, one has 



„ N N-1 

v^p) < z;vt(^-i)^+("-i)(^-i)(i - tr w n^r' n {^^-\f\{dx, 

^^=2 2<i<i<N 7 = 1 



1 



Zn ^(;V-l)= + (a-l){iV-l)^^_^)a-Y^(p^_^) / _ ^^•^N^-2N+iN-l)ia-l) 



Zn-1 

^ _Z^^iN-lfHa~l)iN-l)^l_^^a-ly^^jy^_^y 

Zn-1 

Stirling approximation formula (8) implies that, for e^*^~^~^^ < t < |, 

VR^{IC\V) < Ct, (26) 

for some constant C > 0. 

Now we bound Vhs{^'^) from below. That is, 

N N-l 



VhsP) = Z^t'^''-'^' f n {^^->^JTll{l-t + tX,-tX,fl[dX, 

"^^1 2<i<j<N j=2 j=l 

> ZNt'^'-WnsiV), 
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where the inequahty follows 1 —t+tXi — tXj > tXi — tXj for all j = 2, ■ ■ ■ , TV. Equivalently, 
VR//5(/C*, "D) > t. Together with inequality (26), one has 



for c = -i > 0. This concludes the optimality of the upper bound for a < 1. 

4 the a-volumes of S and WT 

In this section, we provide answers to the questions (1) and (3) for large N in the sense 
of a- volume. These two questions are related to the ones listed on page 9 in [1], that is, 

• Does the volume of the set of separable states go really to zero as the dimension of 
the composite system grows, and how fast? 

• Has the set of separable states really a volume strictly smaller than the volume of 
the set of states with a positive partial transpose? 

Here, we supply solutions to those two questions in the sense of a-volume for large A^. 

The geometry of V is reasonably well- understood, and hence the volumes (e.g., 
Hilbert-Schmidt, Bures, a) of V can be calculated precisely [21, 23, 25, 28]. However, the 
geometry of S and WT are more complicate, and less information for them is available. 
This makes the calculation of volumes (Hilbert-Schmidt, Bures, a) of S and WT quite 
intractable. Thanks to the remarkable works of [12, 13], one knows that VYiHs{^iT>) 
decays to very quick for large N and Y'Rhs{T'VT ,V) is essentially a constant. In 
[14], the author obtained similar results in the sense of Bures volume. Therefore, the 
probability of finding separable quantum states within quantum states is extremely small 
and the Peres- Horodecki PPT criterion as tools to detect separability is imprecise for large 
iV, in the sense of both Hilbert-Schmidt and Bures volumes. 

The following two corollaries give the estimates of a-volume of S. They are direct 
consequences of Theorems 1 and 2, and the estimates for VRhs{S,T>) implicit in [12]. 

Corollary 3 (Large number of small subsystems) Let a > be a (fixed) constant 
independent of the dimension of l-i. For system l-L = (C*^)®", there exist universal 
constants 04, C4 > 0, such that for all D,n >2, 



VRhs{IC\V)>cVR^{}C\V), 



) 



max{l,a} 



( 



(Dnlnn) 



1/2 



) 



min{l,Q} 



_/\/-l/2+a£, 



< VR„(5,P) < C4 



_/y-l/2+ao 



where ao = \ log£)(l + ^) 



1 

2D' 



log^(D + l). 
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Remark. Recall that for T-L = (C^)*^", = Corollary 3 means that, for fixed 
(small) D and fixed a > 0, VRa{S,V) is bounded above by (up to a multiplicative 

/ / \ min{l,o} 

universal constant) ( ) , and hence VRo(iS, V) goes to zero exponentially 

as 77, — 7- oo. This implies that ¥a{S, N, D) =: ^"^j , the a-probability of finding 
separable states within V (on "H = (C^)®*^), goes to zero super-double-exponentially, 
since VRa(5,P) is already the (A^^ — l)-th root of ¥^{3, N, D). Numerical results 
show that, for even (moderately) large N, Fa{S, N, D) is very small. For instance, 
P2(5,256,2) < 2.1 X 10-1595 with n = 8, P2(5,512,2) < 3.1 x iQ-^^mi ^j^j^ n = 9, 
P3(5,512,2) < 1.1 X 10-43412 ^j^j^ ^ ^ g^j^^ P2(5,243,3) < 1.52 x lO-^soi ^j^j^ ^ ^ 5_ 
Similarly, we list some probability for < a < 1, Po.5(iS, 256, 2) < 8.8 x io~479 ^j^j^ 
n = 8, Po.5('5,512,2) < 3.36 x 10^21102 ^-^^-^ n = 9, Po.2(5, 512, 2) < 6.62 x lO^^^^^ ^j^j^ 
n = 9, and Po.5(>S, 243, 3) < 9.5 x lO^^ssi ^j^h n = 5. 

Corollary 4 (Small number of large subsystems) Let a > be a (fixed) constant 
independent of the dimension of Ti. For system H = (C^)®", there exist universal 
computable constants c^,C^ > 0, such that for all D,n > 2, 

max|l,a} / / i M /2 \ min{l,a} 



< VR„(5,I?) < C; 



jYl/2-l/(2n) / - ""V-. ; _ -0 I JYl/2-l/(2n) 

Remark. Corollary 4 shows that, for fixed (small) n, "the order of decay" of VRq,(i5, "D) 
is between and i;(|-f )min{i,a} as D ^ 00. Therefore, P„(5,iV,D) goes to 

zero with order between /}(i/2-n/2)(D2"-i)min{i,a} ^(i/2-n/2)(D2"-i)max{i,a} as D ^ 00. 
Numerical results show that Pq,(iS, A^, -D) is very small for (even moderately) large iV, and 
its "order of decay" relies not only on D but also (heavily) on n. For instance, if a = 2, 
then P2(5,784,28) < 1.35 x 10-^230 ^jth n = 2, P2(5,900,30) < 8.5 x 10-19028 ^j^j^ 
n = 2, P2(5,729,9) < 4.9 x 10-2^218 ^j^j^ n = 3, and P2(>S, 1000, 10) < 1.74 x 10-^7132 
with n = 3. Similarly, if a = 1/2, then Po.5(>S, 900, 30) < 1.12 x 10-8235 ^j^j^ n = 2, and 
Po.5(5, 729, 9) < 6.8 x 10-12538 ^j^^ n = 3. 

For a bipartite system "H = C^^ ® C^^, any state p on "H can be expressed uniquely 

as 

where {ej}^\ and {/q,}^=i are the canonical bases of C-^^ and C-^^ respectively. Define 
the partial transpose T(p) with respect to the first subsystem as 

i,j a,p 
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We use WT (c V) to denote the set of quantum states with PPT, i.e., p G WT if and 
only if T(p) > 0. Thus the Peres- Horodecki PPT criterion is equivalent to 5 C WT. For 
systems (g) C^ (g) and O one has S = WT. That is, the Peres- Horodecki 
PPT criterion is a necessary and sufficient condition for p being a separable quantum 
state if dimension of T-L less than or equal to 6. 

The following corollary is to estimate the a- volume of WT . It is a direct consequence 
of Theorems 1 and 2, and Theorem 4 in [12]. Note the upper bound is trivial because 
WT CV. 

Corollary 5 (a-Volume of WT) Let a > be a (fixed) constant independent of the 
dimension ofH. There exists an absolute computable constant (depending on a) Cq > 0, 
such that, for any bipartite system 'H = C^ cq< VRaiWT, V) <1. 

Remark. An immediate consequence of Corollaries 4 and 5 is that, for = C^®C^ and 
large D, "the order of decay" of VRc,{S,WT) is between and L)-imm{i,a} 

as -D — 7- oo. The upper bound decreases to as i5 — oo, and hence S has really a 
volume strictly small than WT. Moreover, the conditional a-probability of separability 
given positive partial transpose condition is exceedingly small. In other words, the Peres- 
Horodecki PPT criterion is imprecise as tools to detect separability for large A^. We point 
out that our proof yields co{N,a) (e-i/4/4)max{i,a} ^ (^Q^^gg^maxji.a} as ^ oo (see 

also [12]). Any improvements of this constant are very much appreciated. Let us denote 
the conditional a-probability of S given WT as F^{S\WT, N, D) =: for fixed 

a > and A^ = D". Below, we provide some numerical examples to see how small 
of F^{S\WT,N,D) for (moderately) large A^. For instance, Pi.i(5|PPr, 4096, 2) < 

2.5 X 10-2721940 ^j^j^ ^ ^ Pi.i(5|PPr,6561,3) < 1.82 X 10-12248770 ^j^j^ ^ ^ 

Po.9(cS|PPr,6561,3) < 3.74 x lO"™! ^j^h ^ = 8. 

5 Conclusion and Comments 

In the present paper, we compare several important measures (Hilbert-Schmidt, Bures, 
a- measure) on the set of states V onTi. Roughly speaking, for any fixed a > and any 
measurable subset IC C V, the a- volume radii ratio VRq,(/C,'D) is bounded from below 
by VRH5(/C,r')™^^^^'°^ and from above by VRjf5(/C, 1))™"^^'°^ (up to some universal 
constants). We also compare the a- volume radii ratio VRq,(/C, V) with the Bures volume 
radii ratio VRb(/C,'C) and the /3-volume radii ratio VR^(/C,'D) for fixed /3 > 0. In 
particular. Theorem 3 improves the estimates on VRb(/C, ^^) in [14] . 

Employing these estimates to the set of separable states S and to WT (the set 
of states with positive partial transpose), we obtain both upper and lower bounds for 
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VRaiS,V) and VRa,(VVT,V) for fixed a > 0. It is showed tliat VRa(VVT,V) is 
essentially a constant, while VRq,(iS, I^) is of order (at most) A^"'^ for some k > (and 
so is VRq,(5, ■pPT)). This means that a typical (PPT) quantum state is entangled; 
namely, randomly choosing a (PPT) quantum state in P on "H with (even moderately) 
large dimension, this state is entangled with probability (very) close to 1. This may 
be of importance in analysis of quantum algorithms or quantum protocols that rely on 
entanglement. This also shows that, for not-too-small dimensions, the Peres-Horodecki 
PPT criterion is not a precise tool to establish separability. Consequently, more precise 
necessary and/or sufficient conditions for separability are in great demand. 

In applications, one frequently considers T-i as T-L = (C^)*^", often with D = 2 or 
D = 3. The estimates we obtain show that the a-probability of S in V, or of S in VVT 
is very small even for small values of n. For example, the ratio of 2-volumes is less 

than 2.1 x iQ-^^gs if = 2 and n = 8, and similarly < 1.52 x IQ-^soi if D = 3 

and n = 5. Likewise, the ratio of (l.l)-volumes < 2.5 x io~272i94o d = 2 and 

n = 12, and similarly, v^l{vvT) < 1-82 x io-i2248770 a d = 3 and n = 8. 

Concerning question (4), the situation is much less clear. In fact, to answer this 
question, one needs to prove that, for any bipartite system "H = C"^ (8) C^, there exists 
a constant Co < 1, such that VRq(PPT, "P) < Cq. Similar questions can be asked by 
replacing a-volume with the Hilbert-Schmidt volume and the Bures volume. 
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